Temperature in nonequilibrium systems with conserved energy 
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We study a class of nonequilibrium lattice models which describe local redistributions of a globally 
conserved energy. A particular subclass can be solved analytically, allowing to define a temperature 
Tth along the same lines as in the equilibrium microcanonical ensemble. The fluctuation-dissipation 
relation is explicitely found to be linear, but its slope differs from the inverse temperature T^^. 
A numerical renormalization group procedure suggests that, at a coarse-grained level, all models 
behave similarly, leading to a two-parameter description of their macroscopic properties. 
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Understanding the behaviour of nonequilibrium sys- 
tems with a large number of degrees of freedom remains 
one of the major goals for statistical physics. Many at- 
temps have been made to describe such systems in terms 
of a limited set of macroscopic parameters, like in the 
equilibrium case [jj]. In the context of glasses (i.e. materi- 
als with huge relaxation times), an effective temperature 
-first introduced in a phenomenological way ^i2i,|j|- has 
been defined from linear fluctuation-dissipation relations 
(FDR) , on the basis of mean- field spin-glass models p , 
and was shown to satisfy the basic properties needed to 
define a temperature Jo| . Since then, this FDR has been 
tested numerically B , H, 19. Hol Ull Ha Il3l | and experimen- 
tally [ij, 15, 16, filial in many realistic glassy systems. 
Still, no clear consensus concerning the good definition 
of temperature has emerged since the measured FDR are 
not always line 

Apart from glassy materials, other classes of systems 
having a finite relaxation time (granular gases, non- 
hamiltonian spin models, etc.) can be considered as out- 
of-equilibrium in the sense that their dynamics does not 
fulfill detailed balance (DB). Their steady state cannot be 
described in general in the framework of equilibrium sta- 
tistical physics, even though they can share some quan- 
titative properties with equilibrium systems -e.g. criti- 
cal behavior [lUi^l. A few attemps [ill El I2I Ei] 
have been made to define an effective temperature from 
Jaynes' maximum entropy condition 25l| or from FDR, 
but the interpretation of the resulting temperatures re- 
mains to be clarified. 

Since the breaking of DB plays a key role in nonequi- 
librium systems, it is worth distinguishing several fornrs 
of DB. What is often refered to as DB in the literature 
is a canonical version: 
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where M^(/3|a) is the transition rate from state a to state 
/3. This form is just a simple way to enforce canoni- 
cal equilibrium when defining a stochastic model, hence 
its usefulness for numerical simulations. Still, one could 
wish to define stochastic models in a microcanonical sit- 



uation. In this case, the stochastic dynamics should ob- 
vioulsy conserve the energy, and one can also assume a 
microcanonical DB relation iy(/?|Q;) — W{a\(3). Actu- 
ally this form of DB -also called microreversibility- is 
not just a useful recipe but can be given a fundamental 
interpretation in an equilibrium context, as it is asso- 
ciated to the time-reversal symmetry of the underlying 
hamiltonian dynamics. 

Turning to nonequilibrium systems, one expects on 
general grounds that the dynamics breaks the time- 
reversal symmetry due to the presence of fluxes or dissi- 
pation. So it may be of interest to study the most simple 
nonequilibrium stochastic systems which are deflned by 
relaxing only the microreversibility condition, replacing 
it by a more general DB relation W{f3\a)fa = W{a\f3)fp, 
while preserving the energy conservation. 

In this Letter, we study a class C of nonequilibriunr 
lattice models which describe local redistributions of a 
globally conserved energy. A particular subclass Cg, de- 
fined later on, satisfying DB (but not always microre- 
versibility) can be solved analytically, allowing to define 
a temperature Tth along the same lines as in the equi- 
librium microcanonical ensemble. Moreover, we derive 
explicitely the FDR and show that it is linear, with a 
slope different from the inverse temperature T^ , thus 
questioning the relevance of FDR to define temperature 
in non glassy out-of-equilibrium systems. In addition, a 
functional renornralization group procedure inrplemented 
numerically suggests that any model from class C behaves 
macroscopically like a member of the subclass Cs, indi- 
cating that DB is restored on a coarse-grained level. 

Our models are defined as follows. A real vari- 
able — c» < Xi < (X) is attached to each site i of 
a d-dimensional hypercubic lattice with N sites. The 
(pseudo-)energy E — ^ J2i=i ^"i i^ conserved by the dy- 
namics defined by the following stochastic rules. At each 
time step, a link (i,j) is randomly picked up on the lat- 
tice, and a random number q £ [0, 1] is drawn from a 
syirrmetric distribution i/'(g). Then the variables Xi and 



Xj are updated as: 

x'^^±^qixl+x]), x;. = ±^(l-g)(xf+4) (2) 

The sign ± is randomly chosen with equal probability. 
The different models in C are distinguished by 'ilj{q). 
These dynamical rules can be formulated in terms of a 
master equation with transitions rates VF({a:^}|{xi}). It 
is generaly hopeless to find the stationary solution of a 
master equation unless the DB condition is fulfilled. Let 
us define the subclass Cs, consisting of the models for 
which ip{<l) is chosen to be a beta distribution: 

It can be shown that in this case, DB holds, namely: 
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Technical details will be reported elsewhere |26 
Thus microreversibility is recovered only for 
case similar to r] = 1/2 has been studied in [23. The 
steady-state distribution is readily obtained from Eq. Q 
by taking into account the energy conservation and nor- 
malizing the resulting distribution: 
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with Zn{E) = KnE''^-^ and Kn = 2''^r(ry)^/r(ryiV). 
For r] ^ ^, this distribution is clearly non uniform over 
the states of given energy; one can then expect important 
differences with equilibrium systems. 

In the framework of equilibrium microcanonical ensem- 
ble, a well define prescription exists for introducing tem- 
perature. One considers an isolated system (with con- 
stant energy) and divides it into two subsystems. Tem- 
perature is introduced as a thermodynamic parameter 
which takes equal values in both subsystems. If moreover 
the value of this parameter is independent of the choice of 
the partition, the parameter can be said to characterize 
the whole system. 

We consider a partition into subsystems Si and 1S2 
which can exchange energy while keeping Ei + E2 ~ E 
fixed. The energies Ei and E2 are fluctuating, but in the 
limit of large subsystems, the mean value of Ek can be 
identified with its most probable value i?^ . Generalizing 
the equilibrium procedure, the relevant quantity to com- 
pute is then the conditional probability P{Ei \E) that the 
subsystem Si has energy Ei given that the total energy 
is E. Using Eq. 0, one finds: 
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In the usual equilibrium microcanonical ensemble, 
Zn{E) reduces to the phase-space area fijv(-E') of the 
hypersurface with energy E. The most probable energy 
El is found from d\nP/dEi\E* = 0, which gives: 
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This allows to define a temperature T^^ for each subsys- 
tem k through (we set ks — 1): 
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Thus Eq. O implies T/^ = T^f^. It can be checked that 
the common value does not depend on the partition cho- 
sen |2g|, so that this temperature can be safely said to 
characterize the whole system. 

From the expression of Zn{E), one finds Tth = s/ij, 
where e = E/N is the energy density. Also, considering 
subsystem iSi as very small with respect to 52, but still 
macroscopic, one can derive in a similar way the 'canon- 
ical' probability distribution: 
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Another way to introduce a temperature in non- 
equilibrium systems is to consider generalized FDR. This 
approach has received considerable attention since it has 
been given a precise interpretation in the context of 
glasses p- Still, its applicability for non glassy out-of- 
equilibrium system remains to be clarified, and can be 
tested within the present model. To this aim, an ex- 
ternal field h must be introduced to allow for the def- 
inition of a response function. A natural way to in- 
clude an external field is to add to the energy E a 
perturbing term —h^^Xi] one can rewrite the energy 

■^h = 5 Tlii=i{xi ~ ^)^ ^P to an irrelevant additive con- 
stant. Note that E^ has the same form as E in terms 
of the variables Vi = Xi~ h. The dynamics of the v^'s is 
then defined in the same way as for the XiS in the absence 
of field, which is consistent with the equilibrium proce- 
dure. One then recovers for Vi the canonical distribution 
Eq. ®. 

To compute the response function, we assume that the 
system, subjected to a field, is in steady state for t < 0. 
At time i = 0, the field is switched off. The response is 
defined for i > by x(0 = d < N'^ Y^^x, > /dh\h=o. 
From the canonical distribution, the following FDR is 
derived: 
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Although Eq. (|10|l does not lead at first sight to a linear 
relation between x(i) and {xi(t)xi{0)), some simplifica- 



tions actually occur. Indeed, it can be seen that correla- 
tion functions of the form {xi{t)"xi{0)"^) with odd inte- 
gers n and m are all proportional to the 'hopping correla- 
tion function' $(i) = {N-^ YlZi '^^(O), with (/),(i) = 1 if 
Xi (t) = Xi (0) and (j)i (t) = otherwise. More specifically, 
(xi(i)"2:i(0)™) = (x,(0)"+™)$(i) 113. 

As a result, one has {xi{t)xi{0)) — 2e$(t), and 
{xi{t)/xi{0)) = $(i), so that the FDR ^ can be rewrit- 
ten: 
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yielding a fluctuation-dissipation (FD) temperature 
TpD = 2e, which is different from Tth — s/rj (except in 
the case 11 ^ h for which microreversibility is recovered) . 
This leads us to the question: which of the two tem- 
peratures Tth and Tpu is more relevant from a physical 
point of view? Actually, it could be argued that both 
temperatures may be equivalent up to a redefinition of 
the temperature unit: if Tth takes the same value in two 
subsystems, so does Tpo- Still, this conclusion only re- 
mains valid as long as 77 takes the same value throughout 
the system. It is then natural to test a more general dy- 
namics. Interestingly, DB still holds if one introduces on 
each link a different distribution tpij (q) : 
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where r]i can take a different value on each site i; q refers 
to site i and 1 — g to site j as in Eq. 101 ■ Then the micro- 
canonical distribution Pstii^i}) takes the same form as 
in Eq. |SJl, simply replacing JJ^ jxi^''"^ by nj^;^^'''"^. 
Following the same lines as above, P{Ei\E) is easily com- 
puted and leads to equal values of the temperature T^^ 
in both subsystems, with Tth = ^kl {ifjk {{ifjk denotes a 
spatial average of 77^ over subsystem k). On the contrary, 
the FDR formally keeps the same form as previously, 
and the FD temperature remains related to the energy 
density through Tpo = 2e. Choosing {rji} such that 
(^)i 7^ ('7)2, the equality T^ = T^^^ imphes ei ^ £2- The 
equipartition of energy breaks down, which in turn leads 
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Thus the two temperatures Tth and TpD are not equiv- 
alent up to a change of temperature unit, but differ qual- 
itatively since Tpo does not necessarily take the same 
value in two systems which can freely exchange energy. 
Still, if Tth was not known, one could have tried to argue 
that Tpo is the correct teniperature, in a spirit similar to 
the procedure invoked in |a| for glassy systems. Indeed, 
connecting a new site acting as a thermometer to the ex- 
isting system, one may identify its average energy with 
^r, as done also to define a granular temperature |28l| . 
Interestingly, this yields precisely the same result as Tpp,, 
i.e. T — 2e (assuming a uniform 77). That Tth is different 
from TpD in this model thus means that the tempera- 
ture does not reflect only the average energy, but also 
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FIG. 1; ReDormalized distribution tpL{q) for increasing sizes 
L, in dimension d — 1. Full lines correspond to one-parameter 
fits with beta distributions. Inset: parameter tjl from the fit 
vs. L'^ for d — 1 (+) and d = 2 (o); dashed line is the mean- 
field prediction Eq. H3[l . 



the amplitude of the energy fluctuations. For instance, 
as Tth = e/??, a large value of r/ implies a low temperature 
and correspond 

Up to now, we have considered only the subclass Cg 
where ip{q) is a beta distribution, for which a form of DB 
holds. But what happens for more general distributions? 
In particular, one could wonder whether versions of the 
model with beta distributions are in some sense represen- 
tative of the generic behavior of all the models belonging 
to C. If ip{q) is different from a beta law, no DB holds 
|3l| . and there is no clear way to find analytically the 
steady-state distribution Pst{{xi}). Yet, numerical sim- 
ulations show that even for distributions tp{q) far from 
beta laws, the two-point spatial correlation functions still 
vanish in steady state. This is also consistent with calcu- 
lations made in the 'q-model' for static sand piles, which 
exhibits som 

A standard way to coarse-grain the system is to intro- 
duce a renormalization group procedure for ip{q). Since 
the distribution Psti{xi}) is not known, an analytical 
approach seems doomed from the outset. We thus im- 
plement numerically a functional renormalization group 
procedure in the following way: a large system is divided 
into blocks of L^ sites each. Then one runs the dynamics 
and chooses a link in order to redistribute the energy. 
If both sites of the link belong to different blocks, then 
the energies E^ and E^ of the two blocks are computed, 
leading to the renormahzed value qa = El /{El -\- E^). 
The histogram of the values of q^ obtained over the dy- 
namics is recorded, yielding the renormalized distribution 
ipLil)- If the initial distribution ip{q) is a beta distribu- 
tion, this renormalization transformation is exact (see 
below). This procedure has been used for many differ- 
ent initial distributions and all of them lead to the same 
coarsed-grained result. This is illustrated on Fig. 1 for 
L = 2" (data obtained by iterating n times the renor- 



malization with L = 2), starting from the distribution 
^{q) = §|sin(27rg)| (filled circles). Very interestingly, 
tpL {q) converges asymptotically when L increases towards 
beta distributions with exponent 77^ linear in L"* (Fig. 1). 
We have checked that the values of qr are decorrelated. 
Starting with a beta distribution, riL can be computed 
exactly within our renormalization scheme. Indeed, the 
(one-site) distribution of the local energy Si = \x1 is a 
gamma distribution with parameter rj. Since the e^'s are 
independent variables, the block energy Ei, also follows a 
gamma distribution of parameter t^l = r^L"^, leading to a 
beta law with parameter 77^ for q^. The important result 
here is that beta distributions are recovered numerically 
when starting from an arbitrary ^'(9)? even for moderate 
values of L (Fig. 1). Consequently, an effective value of 77 
can be defined for any i/'(g) as rje = tjl/L''- {L^ ^)] Ve is 
the value of 77 which would give the same coarse-grained 
behavior by using a beta ip{q) with parameter yye in the 
basic kinetic rules. Interestingly, 77^. can be computed 
within a mean-field approximation |26l| , yielding 
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where Var(g) = (q^) — (q)^ is the variance of the distribu- 
tion tp{q). This value is in excellent agreement with the 
numerical simulations (Fig. 1). The above results suggest 
an appealing scenario for the description of nonequilib- 
rium systems with a conserved quantity and short-range 
correlations. Breaking the time-reversal symmetry leaves 
considerable freedom to choose the dynamics, but the 
renormalization group procedure shows that the macro- 
scopic behavior can be described by a single parameter 
77e in addition to Tth- This new parameter 77e essentially 
describes the way a globally conserved quantity is dis- 
tributed among the different degrees of freedom. Its value 
is fixed in equilibrium: rjf,q = \ here, but different values 
could be expected in other models -for instance, 77eg = - 

for £' = i Yli \xiY'- Note that the present approach uses 
a renormalization procedure in a context where there is 
a priori no diverging length scale in the system, i.e. not 
close to a critical point J2Q|| . 

In conclusion, we have shown how to define a mean- 
ingful temperature Tth from the conditional energy distri- 
bution of subsystems, within a class of finite-dimensional 
nonequilibrium models with conserved energy. The sta- 
tionary distribution is generally non uniform over the 
states with given energy. The temperature deduced 
from FDR does not coincide with Tth, thus showing that 
FDR are not necessarily the most relevant way to de- 
fine a temperature in out-of-equilibrium (and non glassy) 
steady-state systems. Finally, a numerical renormaliza- 
tion group approach indicates that DB is restored on a 
coarse-grained level even when this property is not sat- 
isfied microscopically. This renormalization procedure 
allows to define a parameter 77e which encodes the de- 
viation from equilibrium. The macroscopic behavior of 



the model is then described by the two parameters Tth 
and 77e, i.e. one more parameter than in equilibrium is 
required. 
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